Spectral gap induced by structural corrugation in armchair graphene nanoribbons 
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We study the effects of tlie structural corrugation or rippling on the electronic properties of un- 
doped armchair graphene nanoribbons (AGNR). First, reanalyzing the single corrugated graphene 
layer we find that the two inequivalent Dirac points (DP), move away one from the other. Other- 
wise, the Fermi velocity vf decrease by increasing rippling. Regarding the AGNRs, whose metallic 
behavior depends on their width, we analyze in particular the case of the zero gap band-structure 
AGNRs. By solving the Dirac equation with the adequate boundary condition we show that due to 
the shifting of the DP a gap opens in the spectra. This gap scale with the square of the rate between 
the high and the wavelength of the deformation. We confirm this prediction by exact numerical 
solution of the finite width rippled AGNR. Moreover, we find that the quantum conductance, cal- 
culated by the non equilibrium Green's function technique vanish when the gap open. The main 
conclusion of our results is that a conductance gap should appear for all undoped corrugated AGNR 
independent of their width. 

PACS numbers; 73.22.Pr, 73.23.Ad, 72.10.Fk 



I. INTRODUCTION 

The spectacular interest that have been raised from 
the recent isolation of an atomic Carbon layer, the 
graphene^, is based on the unusual dynamics the elec- 
trons have in this material. In graphene, the electrons 
behaves as massless relativistic particles giving rise, for 
example, to a sequence of Hall plateau quite different 
than the one observed in 2D electrons systems confined 
in semiconductor heterojuntions. Also the transport 
properties are remarkable, the mobility of electrons in 
suspended graphene could be even higher than in any 
known semiconductor—. Even the existing of graphene 
as a two-dimensional atomic crystal and its stability un- 
der ambient conditions is an surprising fact. Accord- 
ing to the Mermin- Wagner theorem, there is not should 
be long-range crystalline order in two dimensions at fi- 
nite temperature. Even more, fiexible membrane em- 
bedded in three-dimensional space should be crumpled 
because of long- wavelength bending fluctuations. How- 
ever these fluctuations can be suppressed by anharmonic 
coupling between bending and stretching modes. As a 
result, single- crystalline membranes can exist but should 
be rippled^. Indeed, ripples were observed in graphene^. 
It has been proposed that they should play an important 
role in its electronic properties'^'^. In particular, intrinsic 
rippling has been proposed as one of the possible mech- 
anism for electron scattering to explain the variation of 
the resistivity with the number of charge carrier experi- 
mentally seen in graphene^. 

Between the research areas of graphene, a very impor- 
tant one is the study of nanoribbons where the sheets 
are catted with a particular pattern to manage the elec- 
trical properties. Depending on the type of the border 
edges they can be either in, namely, zigzag (ZGNR) or 
armchair configurations. It is known that a tight bind- 
ing model predict that ZGNRs are always metallic while 



AGNRs can be either metallic or semiconductor, depend- 
ing on their width*. However, experiments suggest that 
AGNR are always insulators with a gap that scale with 
the inverse of its width^. The effect of the corrugation 
on the electronic properties of graphene has been studied 
in previous works by a tight binding type model^*^ and 
by ab-initio LDA calculation^^. The rippling induce a 
nonuniform gauge field whose effects has been analyzed 
by analogy with an applied magnetic field. In this sens 
a pseudo Landau levels (LL) were predicted. Regarding 
nanoribbons the effects of the rippling on the conduc- 
tivity were analyzed in a model which include also the 
effect of the disorder produced by charged impurities 
Although this approach is quite realistic it can not iso- 
late only the effect of the rippling in order to known the 
contribution of each perturbation separately. Therefore, 
in this work, our mainly purpose is to analyze the ef- 
fect of the corrugation on the electronic structure and 
the transport properties of AGNRs. First, as a necessary 
previous step we study the electronic properties of cor- 
rugated graphene layers. We reanalyze the condition for 
the appearance of a flat band associated with the zero 
LL. We show that an strong rippling on the sheet is nec- 
essary to produce such a flat band. Then, we focus on 
the armchair border type nanoribbons. Particularly we 
analyze the case of the zero gap band-structure AGNRs. 
We show that due to rippling an spectral gap is open 
in otherwise conducting ribbons. We subsequent analyze 
the quantum conductance by the non equilibrium green's 
function technique^'^ (NEGF) and show that the opening 
of the gap manifest in a insulating behavior of the un- 
doped samples. 

The paper is organized as follows. In Section |ll] we 
provide a detailed explanation of the corrugation model 
adopted and we study the effect of the corrugation on the 
electronic spectra of a graphene sheet. Then, in Section 
mil we study the effect of the rippling on the electronic 
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spectra and the conduction properties of metallic AGNR. 
Finally, in Section ITVl we present the conclusions of our 
paper and we discuss it implicants. 



II. THE EFFECT OF THE RIPPLING: THE 
GRAPHENE SHEET 

Let us start by studying the effects of the corruga- 
tion on the electronic states for the case of 2D graphene 
layer. This question was examined in a series of previ- 
ous papersiSiii. As it will be the starting point of our 
study of a corrugated nanoribbons in the next Section, 
we reanalyze this problem in the present Section. This 
will be presented after a detailed description of the model 
adopted. 



bor (Aa) we have: 

i = to H — -aH^a 
a 

Aa = y^a + {z{x) - z{x')y - a (2) 

where to ~ 2.66eV is the hopping parameter of the un- 
deformed graphene which is taken as a unit of energy in 
the present work, a = ~ 2, whose value is taken 

from Ref. The rippling also produce a bending of 

the p-orbitals. However it has been shown that the rela- 
tive change of the hopping due to the bending is weaker 
than the one due to the change of the bond length^^ . For 
simplicity, we neglect this effect in the present paper. 

The tight binding Hamiltonian in the presence of the 
corrugation becomes: 



A. The model 

We describe the electronic properties of graphene layer 
by means of a tight-binding model with one 7r-orbital for 
each Carbon atom and nearest-neighbor hopping between 
them. Carbon atoms arrange in a honeycomb lattice. It 
is not a Bravais lattice but can be constructed from the 
hexagonal lattice by putting a basis of two atoms. The 
hexagonal lattice could be constructed as R' = I^sli + 
lySL2 with Ix and ly integers, ai = |(3,-\/3) and a2 = 
I (3, — -/S) are the primitive vectors. The basis is given 
by the atom A sited at R' and B at R' + f (1, ^3). a 
^ 1.41 A is the lattice constant and hereafter will be taken 
as unit of distances. 

The corrugation was included as a sinusoidal function 
which module the z coordinate of each lattice site. The 
X and y coordinates remain into its values in the hon- 
eycomb lattice. For simplicity, we assume that rippling 
only depend on the x direction (see fig. [T]for the case of 
a nanoribbon). Then, the adopted z{x) function reads: 

zix) = hosini — a; I (1) 

where /iq is the amplitude and r the period of the rip- 
pling. This simple model for the corrugation has been 
used in previous worksi2*ii. As a superposition of func- 
tions of the type ([T]) with different wave vectors and dif- 
ferent amplitudes could generate a quite general corruga- 
tion, our study should be taken as a starting point for the 
effect of a general corrugated function. Evermore as we 
point to the study of nanoribbons which by definition are 
much longer than wider the one-dimensional character of 
the rippling function is justified. 

The modulated structure induce an spacial variation 
of the hopping parameter t, being now dependent on the 
distance between the atoms. Developing t up to first or- 
der in the perturbed distance between the nearest neigh- 



H = t[zBi+s - ZAi\{c\A^^c,^,B.a + h-c) (3) 

l,S,(j 

where I ~ {lx,ly) denotes the integer coordi- 
nates in the hexagonal Bravais lattice and 5 = 
{(0,0), (-1,0), (-1,1)}. 
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FIG. 1: (a) Schematic view of an armchair graphene nanorib- 
bon. Vertical lines separates the regions adopted for the con- 
ductance calculation described in Sec lIII Bl trough the Hamil- 
tonians displayed on top. The central region possesses ria; = 4 
and Uy = 3, listed as indicated, (b) Side view in the x — z 
plane of the displacement of Carbon atoms due to rippling, 
is the amplitude and r the period, a is the lattice constant of 
graphene. Vertical dotted lines from (a) to (b) gives reference 
for the actual positions of Carbon atoms in the corrugated 
nanoribbon. 
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B. Dirac equation 

A low energy Hamiltonian could be obtained from ([3]) 
expanding the Fourier transform of the electron opera- 
tors around two of the inequivalent points where the dis- 
persion relation vanish, namely the Dirac points (DP). 
We choose these points as given by K ~ and 

K' = if (1) — "^)- In addition we assume a smooth vari- 
ation of z{x\ i.e. r >> a. The low energy physics 
is described by the sum of two Dirac Hamiltonians for 
massless particles in a pseudo-magnetic field. It is given 
by: 

H = VF J dr^{^\{r)(T* ■ {-iV + A{r))^i{r) 

+ ^{r)(T ■ {-iV - A{r))^2{r)) (4) 

where vp — is the Fermi velocity, cr = {ax,<Jy) is 
the vector of the Pauli matrices and cr* its complex con- 
jugate. The two component field is = (c^AiCis) with 
Ci,A{B) creation field operator over the A(B) sublattice 
The index i = 1, 2 refers to states with momentum near 
K, K'. The gauge vector A{r) = {Ax,Ay) were induced 
by the rippling and is given by: 



. = 

y - = ~|(/iog)'(cos(^x) -i-i) 

= Ay + ko (5) 



ko = —^{hoq)^ is the k = Fourier component of Ay and 
g — ^ the wave vector of the rippling. Note that as the 
time reversal symmetry is not broken by the rippling, this 
implies that the coupling of with A{r) have different 
sign than the one of '3/2. This fact will be important 
in the study of the nanoribbons undertaken in the next 
Section. 

The Dirac equation in a magnetic field corresponding 
to the one particle problem generated by eq. @ at i — 1 
could be exactly solved for the zero energy eigenstate. By 
following similar steps as in Ref. [ij we find the following 
two degenerated eigenfunctions: 



«'+o(a;,y) 



*-o(a;,y) 



N 



yLyLx 
N 



(6) 



with F{x) = ^hQqfiin{2qx) and L^^y the system length 
in both direction. is a normalization constant given by 

N — [/o(— ^)]~^ with /o the modified Bessel function 
of zero order. 

The low energy eigenstates of (|4]) can be studied by de- 
generated perturbation theory over the zero energy states 
^'+o(a;, y) and ^^{x^y). The potential Ay is periodic 
with period ^. Therefore the solutions of (U) fulfill the 




FIG. 2: Shifting of the DP |fco| normalized to the y coordi- 
nates of this point {Ky) as a function of Iiq/t (solid line). 
We also shown ^ giving the reduction of the Fermi velocity 
(dotted line). 



Bloch theorem. They will have the form: 



(7) 



where n is a band index and k^ belong to the Brillouin 
zone of the imposed periodic lattice i.e. — — <kx < —. 
Otherwise there is not restrictions on the values of ky. 
Note that ^ is of the form of ([7]) with a quasi-momenta 
of the zero energy eigenstates given by fco = (0, — fco). 
and function ipo.k^ix) given by: 



V'o,fe(a;) 



A^ 



^jLyLx 



5oe 



x) ) 



(8) 



flo and 6o satisfying + 6§ = 1. In ([7]) the Bloch 
function fulfill ^n,k^{x + §■) = i^n.k^{x) and satisfy 
H{k)ipn.k = EnMipn.k, with H (k) given by: 



H{k) = Ha{k)+H,nt{k) 
Hoik)=VF(T* ■ (-^V + l(r)) 

-ffmt(fc) = Vf(T* ■ (fc + fco) 



(9) 



As Ho{k)ipo,k{x) = we identify HQ{k) as the unper- 
turbed Hamiltonian. The approximated near-zero energy 
eigenstates could be obtained by diagonalizing Hint in 
the subspace of the degenerated zero energy states. Hint 
proyected onto this subspace is a 2 x 2 matrix called h 
and given by: 



^'^ I fcx - i{ky + fco) 
= vfct* ■ (k + fco) 



(10) 



vf — N'^vf is a renormalized Fermi velocity. The eigen- 
values of PU)) are: 



Ei{k,,ky)^±VF^Jkl + iky + ko)^ (11) 

We see that the main effects of the corrugation are: 
1. DP sited in K have been shifted to K — ko. 
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2. vp decrease by increasing rippling. 

If we had studied the behavior near K' by analyzing the 
second term in (j?]) we would have obtained instead of eq. 
(ITT1) the resuh: 

E2{K,ky) = ±VF^kl + {ky~ko)^ (12) 

Therefore K' moves opposite than K. As fco is a negative 
quantity the two DP away each other in presence of the 
rippling. In Fig. [5] we show the relationship and ^ 

as a function of — 




-0.02 



FIG. 3: Low energy band structure of uncorrugated (blue) 
and corrugated (red) graphene layers. The corrugation pa- 
rameters are t — 72 a and ho equals 0.0, which correspond to 
the Dirac cone, and 6.0 in units of a, respectively. This re- 
sult qualitatively agree with the obtained by means the Dirac 
equation including the pseudo-magnetic field. 

For completeness, let us provide a comparison, just at 
a qualitative level, of the above findings with the band 
structure of corrugated graphene layers. Since we are try- 
ing to detect the displacement of the Dirac points, first 
we must compute the band structure of the uncorrugated 
case. Note here that since the adopted ripple develops in 
the X axis direction with a period equal to r, as can be 
seen in fig[T] (b), then the unit cell must be extended up 
to include the whole rippled region from where the hop- 
ping becomes repeated again. This enlargement of the 
unit cell produces thus a reduction of the k^ components 
of the Dirac points K and K' due to band reflexion at the 
border of the original first Brilloun zone. As can be seen 
in the figlSj where we show the low energy spectrum ob- 
tained for corrugations whose period is t = 72 a and with 
amplitudes ho 0.0 and 6.0 in units of a, the corrugation 
effectively shifts the Dirac point, located originally at K, 
in the ky direction. Also, in addition to this shift and 
in accord with the earlier predictions, the observed small 
slope of the dispersion bands of corrugated graphene im- 
plies a lower Fermi velocity of electrons near the Fermi 
level. It should be mention here that this result is related 
with the appearance of Landau like levels (LL) induced 
by rippling and studied in related workai^. Note that 
from Fig. [5] we see that large values of ^ are required 



to obtain a total flattening of the band. This could be 
interpreted as that the predicted zero LL appears only 
when the rippling is strong enou gh. This conclusion is 
consistent with the results of Ref. Tl. In this work it was 
shown that as a result of the relaxation of the structure, 
the corrugation decrease and the flat band disappears 
from the spectra. 

III. THE EFFECT OF THE RIPPLING: 
ARMCHAIR GRAPHENE NANORIBBONS 

Let us now turn into the main subject of this work 
which is the study of the effects of corrugations on the 
spectra of graphene nanoribbons. As we have mentioned 
in the Introduction, a tight binding based description of 
the electronic properties of the ribbons predict that, de- 
pending on their width and border edge type, they can 
be metallic, with zero band gap, or insulators^. For ex- 
ample, for armchair edge (figH]) there are specific values 
of the width has for which metallic behavior is expected. 
Otherwise they should be insulators. In spite of that, 
the experiments shown that graphene nanoribbons are 
insulators independent of their width^. Though device 
fabrication process does not give atomically precise con- 
trol on the type of edges, this seems to show that a pure 
tight binding model is not enough to describe the GNR 
as it is for the sheets. Both, electronic correlationi^ and 
edge disorder-iSiii, have been proposed as possibles mech- 
anisms that give rise to this behavior. In this Section we 
show that rippling could give an additional source for the 
electronic gap, at least for the armchair edge type. 

A. Dirac equation with boundary conditions 

As for the infinite graphene sheet, it is possible to study 
the low energy electronic structure of graphene nanorib- 
bons by a Dirac like equation, but in this case adequate 
boundary conditions must be imposed^. For armchair 
edges, the valley states near the DPs, K and K', get ad- 
mixture and metallic behaviors are obtained for certain 
widths. We have described in the previous Section that 
the rippling produce a shift of the two inequivalent DPs. 
As we will shown below, to include only this shifting is 
enough to obtain a drastic change in the electronic be- 
havior of AGNRs due to corrugations. 

The wavefunction for sub-lattices A (B) are therefore 
given by: 

$A(s)(r) - e'(^-*=°)"-*i^^(s)(r) 

+ e'^<^'+'=°>'-vl/2^^(5)(r) (13) 

where are the components of the spinor wave 

function for states near K and 5*2, the ones near K' . 
For AGNR the wave function should vanish at y = and 
y = Ly. 
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If we neglect other effect than the shifting of the DP, 
the translational symmetry is preserved and imphes that 
*i(2),A(B)('") could be written as: 



*i(2),A(B)(r-) = 



Vl(2),A(B)(y) 



(14) 



The solutions of the Dirac equation with the previous 
stated boundary conditions (BC) has now the form: 



02, B = -e 
with the energies given by: 



(15) 



(16) 



The values of fc„ fulfills sin[(fc„ + Ky — ko)Ly] — and 
are given by: 



We have fitted the calculated gap using the least squares 
method for a function of the type A = C(^)^ obtaining 
excellent agreement with correlation coefficients near one 
with a precision of 10~^. The obtained fitting is shown 
in Fig. |4]with dotted line. In accord to the stated above 
results, we have observed that by including higher values 
of ho/r the fitting become less precise. 




kn = —-Ky + ko (17) 

Ly 

In absence of kg , fc„ could vanish giving rise to a gapless 
spectra as is seen from Eq. ([T5)) . This is the case when 
the width of the ribbons is Ly = 3{ny — l)\/3a with 
Hy integer. However when ko is present the condition 
for gapless spectra could not be fulfilled in general. Even 
more for the Ly and Uy giving a gapless spectra in absence 
of the rippling we now have a dispersion of the form e = 
ztvp^k'^ + /cq and therefore a gap given by: 

A = 2vFka - (^^^ (18) 

Thus, from the previous consideration we expect a gap 
for AGNR due to rippling that scale as (^) . A more 
accurate treatment of the rippling will correct the con- 
stant multiplying the scaling law and possible produce 
correction to this scaling for large enough rippling val- 
ues. 

To confirm the previous predictions we have calcu- 
lated numerically the eigenstates of corrugated armchair 
graphene nanoribbons. This was performed following the 
steps described at the end of Sec. [Tl] taken into account 
here the armchair ribbon geometry displayed on fig. [T] 
As we were interested in detecting the gap and in ana- 
lyzing its dependence with the corrugation, we have ex- 
plored a wide set of rippling parameters keeping fixed the 
width of the ribbons in the values in which they possesses 
metallic behavior being uncorrugated. Figure |4] shows 
the zero energy gap as a function of ho/r for the various 
sets of parameters displayed on the plot, for a width given 
by Uy = 8. Note that for the small values of /iq/t included 
in this figure, the gap does not depende on /iq and r in- 
dependently but on the ratio between these quantities. 



FIG. 4: Quadratic behaviour of the zero energy gap as a func- 
tion of ho/r in corrugated armchair graphene nanoribbons. 
Values obtained from band structure calculations for the dif- 
ferent corrugation parameters sets displayed on the plot. The 
dotted- line is a fitting to a quadratic law (see text for details). 



B. Quantum Conductance 

Finally, we analyze how all the above findings affect the 
electronic transport properties of the armchair graphene 
nanoribbons. For this purpose, we employed the non- 
equilibrium Green's function formalis m^'^d^ which pro- 
vides an useful method to study the quantum conduc- 
tance in nanoscopic systems. 

Before a discussion of the results let us present a 
briefly description of the method and the implementation 
adopted. We considered one corrugated unit cell, as cen- 
tral region He, connected to two semi-infinite outstand- 
ing graphene nanoribbons flat leads, all having width 
length Uy — 8 which corresponds to a zero band gap 
structure in the uncorrugated case. The hopping param- 
eter in both leads Hamiltonians, and Hji, being flat, 
is to. Hence, the whole system is then described by 

H = Hc + Hr + Hl + hLc + hLR, (19) 

where Hlc and hLB. are the hopping terms from both 
leads to the central corrugated portion of nanoribbon, 
which we set equal to to (see figlUa)). Within this for- 
malism, the Landauer conductance G{E) of the system, 
in a zero bias approximation, is expressed as 

G{E) = ^TiE) (20) 
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where T{E) is the transmission function given by T{E) = 
TT(TL{E)gciE)rRiE)gl,iE)), with r, = -2Jm(S,(£;)) 
{£ ~ L,R) being the couphngs of the corrugated 
graphene nanoribbon to the leads, and Qc{E) = [E — 
He — El — the total Green's function including 

the leads self energies, and S^j. These self energies 
must be calculated through the leads surfaces Green func- 
tions (SGF). Although it was reported a very time effi- 
cient way to compute the SGF of graphene nanoribbons 
leads (See Fig. 2 of Ref. for the sake of simplicity 
we implemented the recursive iteration procedure, given 
by gi ^ {E — Hi — geTg)'^ , where He is the unit cell 
lead Hamiltonian and Ti is the interlayer coupling in the 
semi-infinite lead. For the small width nanoribbons stud- 
ied here we have checked that a fast convergence is ob- 
tained. In the iterative procedure, we set the tolerance 
as 10"^. 

Figures [5] (a), (b) and (c), show the transmission co- 
efficients obtained for armchair border type nanoribbons 
whose period r is 72 a and with amplitudes ho = 0.0, 3.0 
and 6.0 in units of a, respectively. In the left panels of 
the plots we included the corresponding low energy bands 
structure computed numerically following the description 
given at the end of Section [II] suited for the nanoribbon 
case. The band reflexion at the first Brillouin zone can 
clearly be seen in fig. [5] (a) which correspond to the un- 
corrugated case. Here, and such as is expectable, the 
transmission coefficient is equals to one at the Fermi en- 
ergy. In the fig. [S] (b) the band gap at the Fermi en- 
ergy appears, and accordingly to this the transmission 
becomes reduced. The fact that it do not reach zero is 
produced because the central region possesses only one 
corrugated unit cell. Increasing the amplitude of the rip- 
pling the gap becomes broader as is observed in fig. [5] 
(c) and the transmission goes to zero around the Fermi 
energy. These results are consistent with the obtained 
above by solving the Dirac equation. Here also, a second 
null conductance region appears above a thin perfect con- 
ductance interval at higher energies. In addition to these 
results, obtained keeping fixed and varying ho, we cal- 
culated T{E) for a wide range of corrugation periods and 
amplitudes, finding in all cases the appearance of a gap 
at low energies around Ep. 



IV. CONCLUSIONS 

After having analyzed the Dirac equation includ- 
ing the pseudo-magnetic field induced by the corruga- 
tion, and having explored numerically a wide realis- 
tic range of rippling parameters sets, we have found 
that the corrugations produces a gap in the low en- 
ergy electronic spectrum of otherwise conducting arm- 
chair graphene nanoribbons. We have found that this 
gap scale quadratic with the rate between the high and 
the wavelength of the deformation. Accordingly to this, 
the quantum conductance G{E) of the ribbons, calcu- 
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FIG. 5: Transmission coefficient of corrugated armchair 
graphene nanoribbons. rix = 24, = 8 and ho = Q (a), 
fto = 3 (b) and ho = Q (c), in units of a . In each left panel we 
included the corresponding low energy band structure which 
serve to a better analysis of the behavior of T{E). 



lated by the NEGF formalism, vanish around the Fermi 
energy of undoped corrugated samples. 

These results was preceded by an analysis of the cor- 
rugation effects over the electronic spectra of graphene 
layers. Therein, we have found that the corrugation just 
shift the Dirac points and renormalizes the Fermi veloc- 
ity. These findings agree with the results of numerical 
calculations. 

In view of this, we conclude that the corrugations are 
an important source of gap in the spectra in addition to 
the electronic correlations to which typically is assigned 
this effecti^. Although our study was performed model- 
ing the corrugations by single sinusoidal functions, since 
any lattice deformation could be approximated quite well 
by an adequate sum of sins, we strongly believe that 
the results obtained here are extensible to general cor- 
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rugations. Moreover, our findings could be revealed in Acknowledgments 
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